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We study a plasma instability in semiconductor superlattices irradiated by a monochromatic, 
pure AC electric field. The instability leads to sustained oscillations at a frequency 0J2 that is either 
incommensurate to the drive, or frequency-locked to it, W2 = (p/q) 'W. A spontaneously generated DC 
bias is found when either p or q in the locking ratio are even integers. Frequency locked regions form 
Arnol'd tongues in parameter space and the ratio UJ2/0J exhibits a Devil's staircase. A transition to 
chaotic motion is observed as resonances overlap. 



I. INTRODUCTION 



Miniband transport in semiconductor superlattices 
(ssl) is a fertile ground for observing various nonlinear 
transport phenomena [l|. Miniband SSLs subject to in- 
tense, monochromatic THz-radiation exhibits a variety of 
such effects including dissipative chaos [2n3 , and genera- 
tion of a quantized spontaneous DC bias [5(. A sufficient, 
although not necessary, requirement for the appearance 
of such novel dynamics is the presence of instabilities in 
simpler types of motion. A well-studied example of this 
is the instability occurring near the conditions of neg- 
ative differential conductance (ndc). This is known to 
lead to the generation of a nearly quantized spontaneous 
DC bias via dynamical breaking of symmetry but 
also domains of different electric field strength that in- 
validate the used models [l(| UK ■ The conversion of pure 
THz AC input to DC output, or rectification, has immedi- 
ate applications to detection of the radiation, and there- 
fore finding instabilities that might lead to such broken 
symmetry is of great interest. 



In Ref. [rj an instability different from that observed 
at ndc was reported. A weak probe AC field with a fre- 
quency that is an irrational multiple of the pump fre- 
quency was found to be divergent provided the degree 
of nonlinearity was sufficiently high. Notably this insta- 
bility occurs outside the regions of NDC implying the ab- 
sence of domains. In this paper we report our findings on 
the same phenomenon but using a more elaborate model 
that allows us to account for all harmonics in the pres- 
ence of nonlinearity, study the problem in terms of the 
applied field amplitude Eq and frequency u, and impor- 
tantly consider the complex dynamics that arise following 
this instability. Our goal is twofold, firstly to understand 
the dynamics of miniband electrons with strong nonlin- 
earity and secondly, apply the instability for the detection 
of THz radiation via, for instance, rectification. 
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The variables V and W are the average velocity and ki- 
netic energy along the superlattice axis of the electron en- 
semble following a distribution given by the Boltzmann 
transport equation. E is the total electric field inside 
the superlattice and includes the externally applied sinu- 
soidal and the self-consistent fields, 



E(t) = E cos ut + E sc (t). 



(2) 



The first two equations of Eqs. (JTJ arc the well-known su- 
perlattice balance equations for a sinusoidal miniband[l|, 
1 1 31 ] with a miniband width of Ao, superlattice period a, 
and equilibrium average energy W cq , and with dissipation 
modeled using constant phenomenological relaxation rate 
7 = 1/t. Constants e and £0 are the electron charge and 
relative permittivity, respectively. The model includes 
the effects of displacement currents 0, [HI via Eqs. @ 
and (|Tc|) . E sc connects the total electric field to the elec- 
tric current, and thus introduces an additional degree of 
nonlinearity to the problem. The strength of the current 
to E'-field coupling is proportional to the electron den- 
sity N and the maximum velocity allowed by the super- 
lattice miniband, Vmax = aAo/2h. This is conveniently 
expressed by a single parameter, the plasma frequency 
w p i: 
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The product ui v \t then controls the balance of dissipation 
and nonlinearity. 



II. PLASMA INSTABILITY 



As our model, we will use the sinusoidal miniband su- 
perlattice balance equations @, [H, [TH, [l4| with the self- 



Our focus are the instabilities occurring in the symmet- 
ric limit-cycles of Eqs. (JTJ. By symmetry we here refer to 
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Figure 1. [Color online] Types of fixed-points of the map G, 
(a) WpiT = 3, (b) Wpir = 12. For all cj p it the NDC-instability 
is present, visible here as the regions where the fixed-point 
is a saddle [gray]. Equilibrium is a stable focus-node [white] 
for the rest of the shown parameter range, except near the 
borders of NDC-instability, where it is a node [green]. In (a) 
the NDC type is the only instability, but in (b) also the focus 
can destabilize, resulting in the region where the point is an 
unstable focus-saddle [blue]. 



the invariance of the equations under the transformation 
S: 

S : (t, V, W, E sc ) ^{t + T/2, -V, W, -E sc ). (4) 

A symmetric limit-cycle is then to be understood as a cy- 
cle that is invariant in S: (V(t+T/2),W(t+T/2),E sc (t+ 
T/2)) = (-V(t),W(t),-E sc (t)). Symmetric cycles rep- 



resent a basic type of dynamics that do not support ef- 
fects such as spontaneous generation of DC bias. Our 
interest in them lies in the fact that a loss of stability 
in such a cycle would imply a transition to some differ- 
ent types of dynamics, provided no other symmetric cy- 
cles are stable for the same parameters. The stability of 
symmetric cycles can be studied by considering the fixed- 
points of the map G: G(X(t)) = SX(t + T/2), where 
X(t) = (V(t),W(t),E sc (t)), and S = diag(-l, 1, -1). 
Clearly, a fixed point of G corresponds to a symmet- 
ric cycle of Eqs. (P). To map the parameter-space re- 
gions where a symmetric limit-cycle loses its stability, we 
have computed the fixed-points (equilibria) of G which 
we then characterize according the eigenvalues of the Ja- 
cobian Jy = dGi/dXj. 

Our findings are shown in Fig. [1] where we have plotted 
the type of the equilibrium on (o>, _E )-parameter plane 
for two different degrees of nonlinearity, lo v \t = 3 and 
^piT = 12. We have also set W eq = — Ao/2. Two types 
of stable equilibria are present for both cases. Typically, 
the equilibrium is a focus-node, implying that the steady 
state is reached via damped oscillations. Additionally, 
also node-type equilibria exist. For these, trajectories fol- 
low exponentially converging, non-oscillatory orbits onto 
the limiting motion. 

Turning to the unstable fixed points, for plasma fre- 
quencies UJ P \T < 8 the only unstable fixed-point is a sad- 
dle. These appear as the exponential convergence to 
a node changes to exponential divergence, and can be 
shown to coincide with NDC. This is shown explicitly in 
App.lAl Since here the NDC appears at zero DC voltage, a 
region of absolute negative conductivity (anc) is always 
associated with its appearance. This leads to dynami- 
cal rectificatio n fl3 . [ill , an effect that is studied in detail 
elsewhere 0, SEE, 17 1. We refer to the saddle- type insta- 
bility as the NDC-instability. This instability persists for 
all Wpir, and occurs approximately for (uj, Ep) such that 
eaE /(huj) is a root of Bcssel Jo function[13|. 

As the plasma frequency is increased also the sta- 
ble focus-node loses its stability, turning into an un- 
stable focus-node in a Hopf-bifurcation. At the Hopf- 
bifurcation point the damped oscillations around a focus- 
node become diverging. This can be identified as the 
plasma instability of Ref. [ID as it describes just such dy- 
namics. Importantly, the it indeed appears outside the 
regions of NDC, therefore suggesting that at least electric 
domains will not form. The lowest degree of nonlinearity 
that is required is lu p \t ~ 8, although values in the excess 
of ~ 10 are required to observe the effect in a reasonably 
large range of parameters. We note that the signature 
of plasma oscillations is visible for a much wider range 
of uj p \T as the damped oscillations corresponding to a 
focus-node -type equilibrium. 

In terms of superlattice parameters, the requirement of 
high WpiT can be achieved for instance for a superlattice 
with n ~ 10 18 enr 3 and Ao ~ 100 meV, with other pa- 
rameters being a = 6 nm, 7 ~ 4 THz (r = 250 fs). These 
yield (oj p1 t)/2tt ~ 8 THz, with values of E ~ 18 kV/cm 
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and ui/2w ~ 4 THz (hjJ ~ 31 meV) required to reach 
the plasma instability region. These values are beyond 
what is reported in the literature for wide miniband and 
highly doped or low 7 super lattices [lH- On the other 
hand, previous theoretical studies using self-consistently 
calculated energy dependent scattering rates have had 
u) p \T well in excess of 10 near miniband centerQ. Due 
to field screening effects, very high uj p \ requires corre- 
spondingly high pump frequencies. For this reason sin- 
gle miniband transport model can become questionable, 
since it requires that hui <C A g , where A g is the gap 
between first and second minibands. 

To gain better understanding of the dynamics and 
physics of the instability, we transform our equations into 
nearly Hamiltonian variables, i.e. into a form that differs 
from a Hamiltonian by a small parameter. We rewrite 
our equations in terms of new variables q, p, and A: 
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Physically, p is the mean value of the electron ensem- 
ble (quasi) momentum while A describes an effective in- 
stantaneous miniband width, reduced from Ao due to 
the spreading of the electron distribution in momentum 
space. More precisely, the fraction A/Ao = 1 — V, where 
V is the variance of the momentum distribution (using 
a definition appropriate for periodic distribution func- 
tions). By the virtue of being the opposite of the vari- 
ance, A/A describes the coherence of the ensemble mo- 
tion, i.e. the concentration (bunching) of the electrons to 
the vicinity of the mean momentum p. The interpretation 
of these variables is rigorously justified in Appendix [Bl 
Substituting the above to Eqs. ([T]) we obtain 
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The new governing equations can be thought of as ef- 
fective semiclassical equations of motion for an electron 
ensemble. A notable difference is that A is here a dy- 
namic variable rather than a constant, and that the self- 
consistent field, acting essentially as a harmonic return- 
ing force, is present via the variable q. 

The Hamiltonian limit, 7 — > 0, of Eq. © has a well- 
studied classical analog. Defining the phase = ap/H one 
sees that 9 follows the equation of a driven pendulum: 
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Previous works have also considered the pendulum limit 
of the super lattice balance equations [H, Q. These have 
included damping via a term — j9 on the right-hand 



side to model the effect of dissipation, while supposing 
A = constant. However, such a system cannot capture 
all the dynamics observed in the third order differential 
equations ©. This is because for the damped second 
order pendulum all phase space areas contract along the 
flow, making Hopf-bifurcation impossible. This implies 
that the interaction between the high- and low frequency 
parts, (q,p), and A, respectively, generates the quasipcri- 
odic oscillations. 

Using Eqs. © we have computed analytically ap- 
proximate conditions for the appearance of the Hopf- 
bifurcation in the weak drive limit. We defer the details 
of the analysis to Appendix[X] Our findings show that the 
Hopf-bifurcation occurs for E > E cr i t , where E„- lt ~ > 
as 7 — > 0, and for frequency 



o R + e>( 7 2 ) < w < o H + 0(7 2 ), 



(8) 



where Or is the amplitude dependent natural oscillation 
frequency of the system, Oh > Or is an upper limiting 
frequency for the Hopf-bifurcation, and C(7 2 ) represent 
terms proportional r ) 2 . Both Or and Oh tend to uj p \ as 
Eq — > 0. The plasma instability is in other words first 
triggered at a small positive detuning from the resonance 
u) = £1r,, and exists for a modest range of parameters 
extending to higher frequencies from resonance, and can 
be brought about by even a small Eq provided that uj v \t is 
large enough. The frequency of the modulation of A and 
(j>, q) oscillation, denoted f2, is small, vanishing at the 
resonance u> = Or and increasing as parameters varied 
away from it. 

Turning to physics point of view, the analytic results 
show that the most significant contributing factor to the 
appearance of the instability is the coupling between the 
slow degree of freedom A to the resonance frequency Or, 
and in turn, the down mixing of the fast oscillations to 
the motion of A. This is in direct analogy to what is 
found for coupled high- and low- frequency oscillators [l9j. 
where a very feedback coupling between slow and fast 
parts was found to lead to a Hopf-bifurcation followed 
by transition to chaotic motion. Physically, A affects 
the resonance frequency via the coherence A/Aq. This 
is easily seen by noting that the resonance frequency is 
prop ortiona l to the frequency of free linear oscillations, 
Wpi^/A/Ao [cf. Eq. 0], essentially a coherence-dressed 
plasma frequency. Thus, the slow motion of A couples 
to the fast oscillations by modulating the resonance fre- 
quency, and thereby affecting a change in the amplitude 
of the p, q oscillations, with strongest response occuring 
near u ~ Or. 

The reverse coupling, i.e. down mixing of fast to slow 
motion, is due to scattering induced loss of coherence, 
described by Eq. (|6a|) . The rate of decoherence increses 
as the mass of the electron distribution is offset from the 
band bottom, since scattering events in the present model 
return electrons to the thermal distribution that is cen- 
tered to p — 0. Given the change in A is essentially adia- 
batic at the present range of parameters, the decoherence 
per drive cycle is then dependent on the amplitude of p, q 
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Figure 2. [Color online] Regions of different types of dynamics, lj p \t = 20. Quasiperiodic motion [blue], frequency-locking [light 
gray and black], or chaos [red] appears as a result of the plasma instability. For tongues with locking ratios p/q, indicated 
by the labels for select resonances, such that either p or q is even [light gray regions] a spontaneous DC- voltage is present. 
Transition to chaos is associated with overlapping of the resonance tongues. The region of period- 1 solutions with rectification 
near eaEo ~ 2Ahw extending to high-frequencies outside the shown range is born out of the first NDC-instability region. 

oscillations, with large amplitudes contributing a larger 
loss of coherence. This way a slow modulation of the fast 
oscillation amplitude is converted back to a variation in 
the coherence. 



III. SYNCHRONIZATION, RECTIFICATION, 
AND CHAOS 

The Hopf-bifurcation introduces a new frequency into 
the system by essentially generating spontaneous oscilla- 
tions that appear as a parametric self-excitation, as can 
be seen from Eq. ([7]) and noting that A/Ao is now a 
slowly oscillating function. It is expected, then, that this 
generates new dynamics, and in this section we discuss 
the various limiting dynamics that we have found in nu- 
merical simulations of Eqs. ([T]). 

We have observed three different types of motion 
as a trajectory evolves away from the unstable focus: 
quasiperiodic oscillations, frequency locking, and chaos. 
In Fig. [2] we have plotted the type of attractor corre- 
sponding to fixed initial conditions of V(0) = 0, W(0) = 
0, i? S c(0) = 0. In this example, trajectories ejected from 
the unstable focus-saddle typically converge to quasiperi- 
odic oscillations, making the attractor a torus. This in- 
troduces a new frequency into the system that is incom- 
mensurate to the drive frequency uj. We denote this fre- 
quency by u)2, and compute it by finding the strongest 



peak in Fourier transforms of E sc that is not an inte- 
ger multiple of uj. The frequency u>2 is related to the 
slow modulation frequency of the previous section by 
U2 — ui ± Q. This relation is exact only for small plasma 
oscillations, and typically it is the minus sign that ap- 
plies (frequency approximately u + ft is present in the 
frequency spectrum, but is much weaker). 

Of particular interest is the application of the 
quasiperiodic oscillations to detection of THz-radiation. 
To this end, spontaneous DC voltages such as those ap- 
pearing at the NDC-instability would be desirable. How- 
ever, the Fourier transforms of the net electric field E 
on the tori show no zero harmonic component. The 
Hopf-bifurcation is super-critical and so the amplitude 
of the plasma oscillations is small near the bifurcations, 
but grows rapidly away from the critical curve. The fre- 
quency uj2 varies between uj p \ and 0, with the frequency 
decreasing as the NDC-regions are approached. 

For decreasing uj the plasma oscillation amplitude in- 
creases. This increases the coupling of uj and uj2 modes 
resulting in synchronization of the plasma oscillations to 
the drive. When in synchrony, the relationship 

puj = qu>2, (9) 

where p and q are integers, holds for a finite range of pa- 
rameters. For oj v \t > 10, the frequency locking regions 
appear as Arnol'd tongues in and around the quasiperi- 
odic regime [20] . On the tongues the period of the motion 
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Figure 3. Main figure: Frequency of the plasma oscilla- 
tions as a function of applied field amplitude Eo, ui p \t = 20, 
lot = 14.1. Black and gray dots correspond to stable and 
unstable phase-locked solutions (lu2/lu rational), respectively, 
while blue points mark irrational frequency ratios. Chaos 
[red] is seen to occur where the curve becomes dense with 
plateaus, e.g. near and between 7:9 and 13:17 resonances. 
Inset: Generated spontaneous DC-field, {E ac ) over the same 
range of Eq. Strong resonances at iU2 = 3/4 -w and LJ2 = 4/5 -w 
are clearly visible as relatively large peaks in the DC bias. 



starts as q times the period of the pump held T, and in- 
creases via period doubling bifurcations as chaotic regime 
is approached. A more detailed picture is provided in 
Fig. O where we have plotted CJ2/W over a range crossing 
multiple Arnol'd tongues. The curve bears the shape of 
a classic Devil's staircase where the plateaus indicate fre- 
quency locking. Significantly, if either q or p is even, or if 
a period doubling has occurred, we observe a spontaneous 
generation of a weak, unquantized DC bias. In the inset of 
Fig. [3] we have plotted the value of DC-electric field across 
the superlatticc, (E) = {E sc }. Non-zero values in the DC 
component can be seen at even valued resonances, how- 
ever, no quantization is present. The fact that no ndc 
is observed, implies that synchronization can be used in 
the detection of intense THz-radiation. We note that 
tongues with locking ratios such that p,(j > 1 are seen 
only emerging inside the region of plasma instability. On 
the other hand, for locking ratios 1:2 and 1:3 the corre- 
sponding tongue survives outside the plasma instability 
regime. Although the initial conditions used in Fig. [^pre- 
fer the period-1 solutions, these synchronization regions 
persist to higher lo than the quasiperiodic oscillations. 
This means that the 1:2 resonance can in principle be 
used for THz detection via AC to DC bias conversion in a 
wider range of parameters. 

The resonance tongues form a border between 



quasiperiodic and chaotic oscillations. For decreasing uj, 
the mode-locked plateaus widen and more of them ap- 
pear, viz. Fig. [2] Accordingly, the staircase starts to 
show distinct fractal character as steps form on shorter 
and shorter length scales. Transition to chaos is then 
observed as nearby plateaus overlap. This is seen in 
Fig. [3] for example near 3:4 resonance and between 7:9 
and 13:17 resonances. Note that the resonance plateaus 
are still present, but unstable and embedded in the 
chaotic attractor. Onset of chaos in this system can 
therefore be seen as following the overlapping of nearby 
resonances [2lj| . implying the collision of the heteroclinic 
cycles connecting the saddle points associated with the 
g-periodic orbit. 

The route to chaos as observed here is a "universal" 
type of route to chaos also in dissipativc, drive oscilla- 
tory systems [22l - [24j ]. and dissipative systems with two 
competing frequencies. These include systems described 
by the AC+DC driven and damped pendulum equation, 
such as Stewart-McCumber [25[ model of Josephson junc- 
tions, and systems reducing to circle maps or similar dis- 
sipative standard-like maps. A notable difference here is 
that the second frequency is self-generated via the Hopf- 
bifurcation, and not e.g. an additional external drive 
term. 



IV. SUMMARY 

We have performed an extensive study of the dynamics 
of the plasma instability in semiconductor supcrlatticcs 
under intense THz-radiation. We found an analog of the 
instability occurring for widc-miniband, highly doped su- 
perlattices that was studied in Ref. [l2|. This instability is 
different from that associated with negative differential 
conductivity, and appears for parameters where electron 
transport is expected to be stable against domains. 

The instability leads to sustained oscillations that are 
quasiperiodic at some u>2, frequency locked whereby L02 is 
fixed to a rational multiple of pump-field, 0J2 = (p/q)tu, 
or chaotic. Importantly from an applications point of 
view, the synchronization of the plasma oscillations to 
the drive, pui = quj^^ with an even p or q, generates a 
spontaneous DC voltage that can be used in detection of 
THz radiation. In particular, the p = 1, q = 2 resonance 
appears particularly strong, and survives for a wide range 
of parameters. 

Transition to chaos was found to occur as resonance 
tongues overlap, showing that the plasma instability is a 
precursor to the chaotic dynamics found in semiconduc- 
tor superlattices in previous works. Indeed, it appears 
that this instability underlies much of the complex dy- 
namics observed in highly-doped semiconductor super- 
lattices. Finally, we also elucidated the physical mecha- 
nism sustaining the plasma oscillations. We found that 
scattering induced down-mixing of high-frequency oscil- 
lations into a slow modulation of resonance frequency of 
free miniband electrons is mainly responsible for the ap- 
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pearance of the instability. 



Appendix A: Details of mathematical analysis 

1. Negative differential conductivity implies 
instability 

Consider a symmetric limit cycle and introduce an in- 
field probe of the form Eq cos(fit) where Q is an arbitrary 
probe frequency and E is the probe amplitude. E is as- 
sumed to be slowly varying. The probe induces a change 
in the current density of the form j j + dj/dE^. So 
being, we can write the time-evolution of E from Eq. (|lc[) 
as 



- - Alt dj ~ 

Eq cos fit — UEq sin fit = —En. 

£o dE n 



(Al) 



Multiplying by cos fit and averaging over T one obtains 
Eq. 



En 
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where oq is the differential conductivity at frequency 
fl, <7q — 2d{j cxp(iflt))/dEn,, > 0, and er = 
d(j)/dE . The damping rate of the probe field is then 
(47r/e )5i{o'n}, and thus, if ctq < 0, the probe is unsta- 
ble, and if Q > 0, the corresponding limit-cycle is clearly 
a unstable focus-saddle. For oo < the local divergence 
is exponential and not oscillatory, making the unstable 
point a saddle. 



2. Plasma instability 

Here we present the details of the analytic formulas for 
the appearance of the plasma instability. For simplicity, 
we consider here the limit eEo <C huj. We approximate 
(p, q) with a single uj harmonic with slowly varying am- 
plitude and phase 



hR . . 
p = sm(w£ 
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We substitute these into Eqs. (|6|), solve for the deriva- 
tives of the new variables, and finally, apply the averag- 
ing method[26[ to get the following equations for A, R, 
and $: 
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R = —TR 
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where Jk is the fcth order Bcsscl function of the first kind, 
A*(R) = 2|W oq |J (i?), T(A,R) = 2 7 |W eq | Ji(R)/RA, 



and flp> is the free, unforced, undamped oscillation fre- 
quency, 



fl R (R,Af 



A 2J 1 (R) 
A~ R ' 



(A7) 



Approximate symmetric limit cycles are given by steady- 
state values of R, $, and A. The phase $ can be elimi- 
nated from the equations, while A = A*(R) and R sat- 
isfies the implicit equation 
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Dynamics near the equilibrium are determined by the 
Jacobian of the left-hand side of Eqs. (|A6|) . Frequency 
of small plasma oscillations is found up to first order in 
7 to be 



fi = 







(i 


fl%\ 
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where the second frequency fls arises from R dependence 
of the SIr, 



n 2 s = tf pl j (R) [MR) -MR)] 



(A10) 



u pl — 2|W eq |Wp 1 /Ao. Frequencies fin and Sis sat- 
isfy Os < ^r for all uj and R, thus, the frequency fl is 
real for uj < fig and uj > J7r. 

Instability regions are given by applying the standard 
Routh-Hurwitz stability criterion to the characteristic 
polynomial of the Jacobian: 



i / n R 

— ID 

4 V w 



Ul 



+ (2 7 -r)(s + r) <o, 

(Alia) 

+ r(s-r)<o, (Aiib) 



where flu > Or and < flR are 7—^0 upper Hopf- 
bifurcation limit and lower nonlinear resonance limit, re- 
spectively, 



fl 



H,M 



n 2 ±2Q 2 pl 



MR) 2 



and 



1 - 



MR) 
RMR)' 



(MR) + RMR)) 



(A12) 



(A13) 



First of the inequalities IA11I is true in the presence of 
the plasma instability, while the second holds for the in- 
stability arising from the nonlinear resonance. In short, 
Eqs. (|A11[) describe two instability bands, one below reso- 
nance (arising from multistability) and second above res- 
onance (plasma instability). Noting that the second term 
in both inequalities is C( 7 2 ), Eq. §E§ follows. 
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Finally, the second term in Eq. (|A12|) . is critical for 
the plasma instability. If it were too small, fin would fall 
below £7r , and the instability would not appear. Looking 
in an earlier stage of the derivation, this term equals 



a 



dA* gOg 
OR dA 



(A14) 



This crucially depends on the strength of the coupling of 
the amplitude of the fast oscillations to A (dA* /dR co- 
efficient), and the coupling A to the free oscillations fre- 
quency (Ott^/dA coefficient). This leads us to conclude 
that the down-mixing of the fast motion in the equation 
for A and the subsequent modulation of the free oscilla- 
tion frequency is the main physical mechanism generating 
and sustaining the plasma oscillations. 



Appendix B: p is the ensemble mean 
quasimomentum, A/Ao coherence 

Here we justify our claim that p as introduced in 
Eqs. (© is the center of mass of the electron distribu- 
tion in (quasi) momentum space, while A/Aq describes 
the coherence of Bloch oscillations, or the concentration 
(bunching) of the electrons to vicinity of the mean mo- 
mentum. We will use the short-hand notations 8 = ap/h 
and A = A/A . 

Let f(p z ,t) be the electron momentum distribution 
function as given by the Boltzmann transport equation, 
where p z is the momentum along the direction of the cur- 
rent, and let overline denote averaging against f(p z ,t), 

(•) = (a/nh) f*% a , a (-)f(p z ,t) dp z . Due to the miniband 
structure, f(p z ,t) is a periodic function of p z , and for 
conveniance we express p z by an angle <j>, <f> = ap z /h, 
— 7r < (f) < 7r. We first note that the naive expectation 
value of the angle <fi, cannot correctly represent the center 
of mass position, and consequently p z is not a suitable 



choice as the mean momentum. This becomes clear by 
considering f(<f>,t) that is condensed to a narrow, sym- 
metric peak at <fi = tt. In such a case, the expectation 
value <f> = 0, in complete opposition to the "true" value 

Of 7T. 

A more appropriate definition for the mean angle, 
which we will denote by </>* , is given via the first trigono- 
metric moment mi of the distribution [2 7| : 



b* = argmi, mi = exp(i<; 



(Bl) 



where argz gives the complex argument of z. Further- 
more, in analogy to the use of variance of distributions 
on the real line, the spread of a circular distribution can 
be characterized by its circular variance V, 



V = 1 - ImJ 



(B2) 



Note that < V < 1. Zero V is equivalent to the entire 
distribution being concentrated to the mean angle (/)*, 
while unity V implies that a well-defined mean angle does 
not exist. 

We show next that <j>* and V coincide with 8 and 1 — A. 
Returning to the definition of V and W, 



A 



V = V ma , x sm(ap z /h), W 



cos(ap z /h). (B3) 



From Eqs. (0 and (|B3[) one finds that 



As'md = siii(ap z /h), Acos8 = cos(ap z /h), (B4) 



or equivalcntly Acxp(i8) = exp(i</>). It then follows 
directly from the definitions, Eqs. (|B1[) and (|B2[) that 
9 — <f>* and V = 1 — A. Scaling angles back to momen- 
tum units, we have p* = tuff /a = H6/a = p, proving our 
assertion that p is the center of mass momentum. 



[1] A. Wacker, Phys. Rep. 357, 1 (2002) 

[2] K. N. Alekseev, G. P. Berman, D. K. Campbell, E. H. 

Cannon, and M. C. Cargo, Phys. Rev. B 54, 10625 (1996) 
[3] J. C. Cao, H. C. Liu, and X. L. Lei, Phys. Rev. B 61, 

5546 (2000) 

[4] Y. A. Romanov and Y. Y. Romanova, J. Exp. Theor. 

Phys. 91, 1033 (2000) 
[5] K. N. Alekseev, E. H. Cannon, J. C. McKinney, F. V. 

Kusmartsev, and D. K. Campbell, Phys. Rev. Lett. 80, 

2669 (1998) 

[6] K. N. Alekseev and F. V. Kusmartsev, Phys. Lett. A 305, 
281 (2002) 

[7] K. N. Alekseev, P. Pietilainen, J. Isohatala, A. A. Zharov, 
and F. V. Kusmartsev, Europhys. Lett. 70, 292 (2005) 

[8] J. Isohatala, K. N. Alekseev, L. T. Kurki, and 
P. Pietilainen, Phys. Rev. E 71, 066206 (2005) 

[9] J. Isohatala and K. N. Alekseev, Chaos 20, 023116 (2010) 
[10] S. A. Ktitorov, G. S. Simin, and V. Y. Sindalovskii, Fiz. 



Tverd. Tela 13, 2230 (1971)Sov. Phys. Solid State 13, 
1872 (1972) 

[11] M. Biittiker and H. Thomas, Phys. Rev. Lett. 38, 78 
(1977) 

[12] Y. A. Romanov, Fiz. Tverd. Tela 21, 877 (1979)Sov. 

Phys. Solid State 21, 513 (1979) 
[13] A. A. Ignatov, E. Schomburg, J. Grenzer, K. F. Renk, 

and E. P. Dodin, Z. Phys. B 98, 187 (1995) 
[14] K. N. Alekseev, E. H. Cannon, J. C. McKinney, F. V. 

Kusmartsev, and D. K. Campbell, Physica D 113, 129 

(1998) 

[15] Y. A. Romanov, J. Y. Romanova, L. G. Mourokh, and 
N. J. M. Horing, J. Appl. Phys. 89, 3835 (2001) 

[16] D. H. Dunlap, V. Kovanis, R. V. Duncan, and J. Sim- 
mons, Phys. Rev. B 48, 7975 (1993) 

[17] J. Isohatala and K. N. Alekseev(2011), unpublished 

[18] E. Schomburg, M. Henini, J. M. Chamberlain, D. P. 
Steenson, S. Brandl, K. Hofbeck, K. F. Renk, and 



8 



W. Wegscheider, Appl. Phys. Lett. 74, 2179 (1999) 
[19] D. V. Shygimaga, D. M. Vavriv, and V. V. Vinogradov, 

IEEE Trans. Circuits Syst. I 45, 1255 (1998) 
[20] Numerical data of Ref. [3 shows traces of regions of 

synchronization. This phenomenon was not investigated, 

however, possibly because the numerical methods were 

not able to detect frequency locking. 
[21] B. V. Chirikov, Phys. Rep 52, 263 (1979) 
[22] M. H. Jensen, P. Bak, and T. Bohr, Phys. Rev. A 30, 

1960 (1984)T. Bohr, P. Bak, and M. H. Jensen, 30, 1970 

(1984) 



[23] C. Knudsen, J. Sturis, and J. S. Thomsen, Phys. Rev. A 
44, 3503 (1991) 

[24] M. Feingold, D. L. Gonzales, O. Piro, and H. Viturro, 
Phys. Rev. A 37, 4060 (1988) 

[25] W. C. Stewart, Appl. Phys. Lett. 12, 277 (1968)D. E. 
McCumber, J. Appl. Phys. 39, 3113 (1968) 

[26] F. Verhulst, Methods and Applications of Singular Per- 
turbations (Springer, 2005) 

[27] K. Mardia and P. Jupp, Directional statistics, Wiley se- 
ries in probability and statistics (Wiley, 2000) 



